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Using wind turbines (WTs) equipped with doubly fed induction generators 
(DFIG) is a popular technology for generating renewable energy. To ensure 
safe operation, prompt maintenance, and better operational reliability, the 
induction generator used in wind energy must be monitored. In this paper, an 
analysis is carried out on stator currents of the DFIG machine in a wind farm 
to identify any blade imbalances in the wind farm. A fault characteristics 
extraction analysis is carried out on the machine stator currents to detect the 
fault in the system. Firstly, the mathematical equation of the DFIG blade 
unbalanced stator current is generated using the DFIG model. Secondly, 
Park's Transformation is used to modify the stator's 3-phase current. Further, 
by evaluating the feature frequency amplitude variation in the squared signal 
by doing a spectral analysis on the stator current vector's squared signal. 
Lastly, a Simulink model for the DFIG is developed. The suggested 
approach analyses the fault signal of the imbalanced blade fault at various 
wind velocities. The outcomes show that the suggested method for 


diagnosing impeller imbalance faults can successfully locate the fault. 
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1. INTRODUCTION 

Early fault diagnosis reduces machine failure by 75%. Timely fault identification is crucial for 
reliable wind turbine (WT) operation. Fault diagnosis tech is critical in power and energy generation by 
making it more dependable. The marine current turbine's rotor blade pitch imbalance defect is proposed in 
[1]. Continuous Morlet wavelet transform and K-nearest neighbour (KNN) machine learning are used to 
classify fault severity in fault diagnosis. MRnet, a novel technique, utilizes photos taken by automated aerial 
vehicles to identify blade defects in WTs. Mask region-based convolutional neural network (R-CNN) and 
MRnet networks are used for fault detection, segmentation, and classification. proposes a methodology for 
frequency control of microgrids using droop control in the generator-side converter (GSC) with battery 
energy storage system (BESS) [2]. The spatiotemporal distribution of natural wind speed and the effects of 
wind shear and tower shadow effect are analyzed along with blade mass imbalance fault in DFIG WTs [3]. 
Frequency spectrum and Lissajous curves are studied and analysed to detect WTs flaws [4]. Zhang et al. [5] 
proposes an image array for recognizing the defect in turbine blades using a neural network. Finally, it 
diagnoses faults in the field-oriented control of doubly fed induction generators (DFIG) using frequency 
spectrum and Lissajous curves [6]. 
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Han et al. [7] combined the unstable streamlined and WT dynamic models to compute a unique 
response based on beam theory in a typhoon for 6 MW and 2 MW WTs. Meanwhile, Rolán and Pedra [8] 
proposed a calculation to investigate DFIG WT. At the same time [9] highlighted the need for a geographic 
information system and presented an improved condition monitoring system, and [10] used a parameter 
estimation technique and an observer approach to diagnose mechanical drive train faults in WTs. 

The use of the artificial neural network (ANN) in WTs is analyzed [11]. The study outlines the main 
challenges and technology gaps. Another study [12], highlights how an inter-turn short circuit (ITSC) fault 
can impact the performance of an induction generator in WTs. The effect of increasing the power converter 
switching frequency of DFIG during fault conditions is investigated [13]. To improve the fault ride-through 
(FRT) capability of DFIG, the switching frequency control parallel resonance fault current limiter (SFC- 
PRFCL) is applied. Additionally, Shohag et al. [14] reviews the various types of damage that WT blades can 
experience and the preventive measures. Lastly, a novel multi-level inverter is developed in [15] to protect 
the DFIG from voltage drops caused by power system issues. 

Zhao et al. [16], discuss WT impeller mass moment computation under five different circumstances. 
Hichem and Tahar [17] presents fuzzy logic for fault monitoring and identification. Yang and Chai [18] is a 
review of fault diagnostics research for wind converters. Qiao and Lu [19] presents a comprehensive review 
of condition monitoring and fault identification methods for offshore WTs. Patel and Patel [20] presents a 
complete study on WT's transient stability. Mishra and Shaikh [21] presents a model for simulating steady- 
state behaviours of a DFIG. 

The researchers [22]-[24] respectively report on rotor unbalance, acoustic emission detection method, 
and PQ theory for permanent magnet synchronous generators. Kusiak and Verma [25] analyses data-mining 
algorithms using data from a large wind farm. An algorithm is introduced that restructures mass and aerodynamic 
imbalances by [26]. Amirat et al. [27] describes a tool for detecting bearing issues and the stator current signal. 

DFIG WTs have different faults such as gearbox faults, blade imbalance faults, and interturn faults. 
The most common fault is the blade imbalance fault. It is necessary to identify these faults early for the 
reliability of the system. Various techniques are applied to identify these faults such as KNN, ANN, CNN, 
and MRnet. The power spectral density method and the order tracking analysis method are proposed for 
detecting imbalance faults by analyzing aerodynamic torque [28]. The spatiotemporal distribution of natural 
wind speed and the effects of wind shear and tower shadow effect are analyzed along with blade mass 
imbalance fault in DFIG WTs [3]. 

In this paper, an analysis is done on stator currents of the DFIG WT to identify any blade imbalance 
conditions in the wind farm. The results of spectrum analysis and the various wavelet decomposition 
components of the signals show the abnormalities of the signal with the blade imbalance fault. Whenever the 
wind speed exceeds its critical limit, the amplitude characteristic frequency Ip’(t) is also exceeded. The defect 
characterization methodology is used to monitor the changes in the magnitude of the characteristic frequency 
I*(t) under the same wind speed condition to identify the severity of impeller imbalance flaws. This analysis is 
carried out with the variable wind speed which changes with respect to the real time. To validate the results, 
wavelet analysis is carried out on the flow of the reactive power of DFIG to determine the disturbances in the 
signal during healthy and faulty conditions. A DFIG Simulink model with blade imbalance is developed in this 
study using stator current characteristic analysis and its coordinate transformation. Additional stator currents are 
evaluated in comparison to their spectrum. The fault characteristic component of the characteristic frequency 
magnitude, Io(t), grows as the wind speed increases. By keeping track of variations in the stator current fault 
characteristic component's characteristic frequency amplitude under the same wind speed conditions, the 
component may accurately assess the sternness of the blade imbalance fault. 


2. DOUBLY FED INDUCTION GENERATOR 

Nowadays modern power electronics converters control the rotor currents in the DFIG to achieve the 
variable speed necessary to maximize energy from the variable wind. Figure | depicts the standard non- 
uniform velocity WT system with a DFIG. Rotor turbine blades, a gearbox, a generator, and a power 
electronics converter are the main components of the DFIG machine. Research by Hammouchi et al. [6] 
GSC, DC link, and rotor-side converter (RSC) are the components that make up the power converter module. 
The common connection is used to couple the stator, GSC, and DFIG to the grid. The dynamic equations of 
DFIG in d-q components applying the Park transformation are given in accordance with the motor sign 
convention and the synchronous reference frame. 


Vsa = (Rs + Lsp)isa + Mpira — Lswsisq = Mwsirq (1) 
Voq = (Rs at LsD)isq a Mpi,g g Lswsisa 7 Mosia (2) 
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Via = (Ry + L-P)ira + Mpisa — LySWsipg — MSwsisq (3) 
Yq = (R, + L,P)irg + Mpisq — LrSwsira — MSQsi sq (4) 
Im = PM(ira isq z i-glsa (5) 


where V,q is stator voltage direct component, V,, is stator voltage quadrature component, V,g is rotor voltage 
direct component, V,g isrotor voltage quadrature component, I,q is stator current direct component, 
I5q is stator current quadrature component, [4 isrotor current direct component, lq isrotor current 
quadrature component, R, is stator winding resistance, R, is rotor winding resistance, L, is stator winding 
reactance, L, is rotor winding reactance, and p is d/dt is the first order derivative operator. 


VAs 


P, e 
<= Rororsde Gidse k 


converter converter 
Figure 1. Considering the motor sign convention, DFIG-based WT with active powers representation 
Im shows electro-magnetic torque, d and q are the direct and quadrature components of the Park 
transformation respectively; and s represents stator and r represents rotor. s=(@s—9 œm)/œs is the DFIG slip. 


Figure 2 shows equivalent circuits that can be used to depict dynamic equations of DFIG. 
Where: þsd=Ls isd, ġsq=Ls isq, þrd=Lr ird, and orq=Lr irq. 


R: -ODs LM L-M "SOba. R Osdsd K-M L-M SOs 
Isd "ltd Isq 


Vsd M Vrd Vsq M Vra 


Direct component Quadrature component 
Figure 2. DFIG equivalent circuit with d-q components 
The instantaneous value of active and reactive power that is switched amongst the grid and the WT 
is given as: 
P = Ps + De = salsa + Vqisa = Vaira + Vegirq 
q = qs + Gr = Vsaisa + Vqisa 


Where the converter's power losses are disregarded, i.e., p-=pr and the GSC is expected to be operate with 
unity power factor, i.e., q,=0. 


3. MATHEMATICAL MODELLING FOR BLADE IMBALANCE FAILURE’S EFFECT ON THE 
STATOR CURRENT OF A WT 
When a WT has three blades and is run at a fixed angular speed of m (rad/s), the blade imbalance 
fault can be said to be present in terms of shaft torque. 


T = To + Ty sin(wmt +) (6) 
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Here t(s) is time; T (Nm) is WT shaft, output torque; To (Nm) is aerodynamic torque; Ty (Nm) is fluctuation 
torque due to blade imbalance fault; and (rad) is early phase angle. 

With its wheels, drive shaft, gearbox, and generator rotor, the DFIG transmission mechanism was 
equivalent to a centred mass symbolised by: 

2Hy 2" =T -Te — Dy@m (7) 
Here Hm is equivalent inertia time constant; T.(N-m) is electromagnetic torque; and Dm, damping 
coefficient~0, at a steady condition, d@m/dt=0. 

In (6) and (7), show the electromagnetic torque oscillates at the same frequency as the shaft torque's 
cyclic oscillation. Consequently, the electromagnetic torque can be explained as follows when the blade is 
unbalanced: 


Te T Teo T Tev SiN(Wmt F ọ) (8) 


Here Teo (Nm) is electromagnetic torque established by the related aerodynamic torque To (Nm) and Tev (Nm) 
is electromagnetic torque formed by the equivalent imbalance torque T, (Nm). 
The isq component in the q-axis direction of the stator current can also be used to express the electromagnetic 
torque Te (Nm): 

3 


Te = F= Yigg (9) 


Where p is No. of pole and s(Wb) is stator winding flux linkages. 
From (10) can be calculated using the DFIG stator flux orientation model: 


aY a 
ae = — Rsisa + Usa 
disd _ LsRy+ LrRs . RrYsd : 
oL, w age bi + age OL, Weisq + Usa (10) 
saa 
~ LL, 


where Sq is stator flux s, d-axis component, isa is stator current is, Usa is Stator voltage us, Rs and R; is stator 
and rotor resistance, respectively, Ls and L, is stator and rotor inductance, respectively, Lm is magnetizing 
inductance, and œs is slip angular velocity. 

To obtain Y,, in the steady state operation of the WT, we need eliminate the u,q in (10). 


. Ls f 
Ysa = Lsisa — o g bsWslsq (11) 


From the stator flux orientation model, the Yq that is the q-axis direction component of the stator 
flux Y, is equal to 0. Thus, 


i Lr ; 
Po = Ysa = Lsisa — OF LsWslsq (12) 


In the steady state case, the stator flux ¥, is constant. The periodic oscillation in Te causes the q-axis 
component of stator current isq to have the same vibration frequency. Thus: 


isq = İsqo + AsigSiN (Wmt + Pq) (13) 


Similarly, when combined with the stator flux Y expression in (12), the d-axis component of stator current 
i,q also shows the same vibration frequency, as (14): 


İsa = isao + Agiqsin (@mt t Pq) (14) 


In (13) and (14), the stator current components isqo and isao are generated by To, respectively. Asia and 
Asiq are the amplitudes of the additional stator currents generated by the imbalanced torque Ty, which are 
proportional to Tev. da and ġq are the initial phase angles of the current isa and isq, respectively. 

The stator current can be transformed from the dq coordinate system to the abc three-phase coordinate 
system when combined with the dq-abc transformation, the calculated a phase stator current is as (15): 
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isa (©) = ip sin(wet + Po) 
+ st Asia cos [Cwe + Wm) t— Pal + Asiq cos [(We z Wm) t— Pal} (15) 


z ~{ Asia cos [(We + Wm) t+ ba] z Asiq cos [Cwe z Wm) t+ Pal} 


where ig = yiio + i240» ho is arctan (isqo /isao ), and We is the synchronous rotation angular velocity. 
Usually, Pa ~ Pq, (15) can be simplified further as (16): 


İsa (t) = lo sin(Wet + Po) + A, cos [Cwe ~ Wm) t— al + Asiq cos [we + Wm) t— Pq (16) 


Where A, = “asd and A, = “aA 

In (16) indicates that the blade imbalance would cause the DFIG stator current to produce a 
harmonic component with frequency fe + fm (we=2 af e, wm=2 afm) (Hz). Although the harmonic 
component can theoretically be utilized to diagnose blade imbalance failure, the actual fault's effect on the 
harmonic component's amplitude is typically minimal. Thus, base frequency current leakage and outside 


noise can readily drown out the amplitude. 


3.1. Method used for doubly fed induction generators system analysis 

Spectrum analysis detects errors by identifying changes in the amplitude of the featured frequency 
of the squared signal stator current vector. To extract fault characteristics, a feature extraction technique 
using Park transform, wavelet and spectrum analysis is intended. Blade imbalance results in harmonic 
components with a characteristic frequency of fe + fm in the stator current. Direct spectral analysis cannot 
extract fault characteristics. Coordinate transformation can analyze stator-phase current fault characteristics 
successfully at consistent and tempestuous wind velocities. Wavelet and spectrum analysis can diagnose 
blade unbalance. See Figure 3 for the steps. 


Acquiring the stator 3 — 
phase current signal i(t) 


' 


Applying Park’s 
Transformation abc - dq 


x 


Envelop Analysis to 
eliminate DC component & 
attain Io (t) 


Using the wavelet to get the 
spectral information and 
extracting amplitude of the 


The amplitude of 
characteristic frequency 
fm at normal operation 


characteristic 


Output: blade imbalance, failure 


Figure 3. Flow chart for fault diagnosis of DFIG WT 


4. RESULT ANALYSIS 

Utilizing Simulink, a DFIG simulation model with per-unit values and a blade imbalance condition 
is created. This model validates the effectiveness of the suggested approach for the stator current fault 
analysis. Park’s transformation methodology is applied to determine the impeller imbalance faults in the 
DFIG. 

The block diagram model is presented in Figure 4. This block diagram consists of a 3-phase source, 
a transmission line of 30 km, and a step-down transformer which is directly coupled with the DFIG WT. The 
dq-axis components as well as the fault characteristic component with fm frequency of the DFIG stator 
current can be observed through the scope under the normal and blade imbalance fault conditions. 
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Figure 4. Simulation model 


4.1. Characteristic analysis of stator-phase current fault 

Direct spectrum analysis on stator current and turbine speed shows that transforming the 3-phase 
stator current from 3-d to 2-d coordinates using Park's transformation can accurately identify fault 
characteristic frequency, even at high wind velocities, and can successfully detect blade: 


ig = Éi Li Li 
s 3a 6 P 6 ¢ 
ma ie i; 
ig = [ol zie 


P(t) =i + ig 


Park's transformation and the squared signal I(t) is used to convert the DFIG stator 3-phase current 
i(t) into two-dimensional space. Park's transformation's main idea is to translate the 3-phase stator current 
into 2-D coordinates. Using envelope analysis to obtain the upper and lower envelopes of I*(t). By removing 
the DC component: 


I (t) = P(t) = aap (t) + 1own(t)]/2 
The frequency of the fault feature amplitude can be determined by doing spectral analysis shown in 


Figure 5. By contrasting the feature frequency amplitude with the variation in the feature frequency 
magnitude during normal action, the savagery of the blade imbalance fault may be assessed. 


Time (seconds) 


Figure 5. Stator current fault characterization analysis 
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4.2. Characteristic analysis of stator-phase current fault 

The characteristics of DFIG under normal condition is shown in the Figure 6. The 3-phase p.u. 
voltage at 575 V bus, 3 phase p.u. current through 575 V bus, and flow of real and reactive power shown in 
Figures 6(a)-(d) respectively. From the obtained results it is analyzed that DFIG operates under blade balance 
conditions. Further, Figures 6(e)-(h) show DC excitation voltage, rotor speed in p.u., and the 3 phase p.u. 
voltage at 25 kV bus, 3 phase p.u. current through 25 kV bus respectively. These results validate the normal 
operating condition of DFIG. 
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Figure 6. The characteristics of DFIG under normal condition; (a) 3-phase p.u. voltage of DFIG at 575 V bus 
under normal condition, (b) 3-phase p.u. current of DFIG at 575 V bus under normal condition, (c) 3-phase 
active power under normal condition, (d) 3-phase reactive power under normal condition, (e) DC voltage of 
DFIG under normal condition, (f) rotor speed under normal condition, (g) 3-phase p.u. voltage of DFIG at 

25 kV bus under normal condition, and (h) 3-phase p.u. current of DFIG at 25 kV bus under normal condition 


To investigate blade imbalance in a DFIG WT, a fault was simulated by changing the pitch angle of 
the blades at 0.4 seconds. The resulting characteristics are shown in Figure 7. A comparison was made 
between the normal and imbalance conditions. Due to blade imbalance condition various changes are 
observed in the results are described in Figure 7. 

The 3-phase p.u. voltage at 575 V bus, 3 phase p.u. current through 575 V bus, and flow of real and 
reactive power shown in Figures 7(a)-(d) respectively. From the obtained results it is analyzed that DFIG 
operates under blade balance conditions. Further, Figures 7(e)-(h) are shows DC excitation voltage, rotor 
speed in p.u., and the 3 phase p.u. voltage at 25 kV bus, 3 phase p.u. current through 25 kV bus respectively. 
These results validate the normal operating condition of DFIG. 
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Figure 7. The characteristics of DFIG under blade imbalace condition; (a) 3-phase p.u. voltage of DFIG at 
575 V bus under fault, (b) 3-phase p.u. current of DFIG at 575 V bus under fault, (c) 3-phase active power 
under fault, (d) 3-phase reactive power under fault, (e) DC voltage of DFIG under fault, (f) rotor speed under 
fault, (g) 3-phase p.u. voltage of DFIG at 25 kV bus under fault, and (h) 3-phase p.u. current of DFIG at 
25 kV bus under fault 


4.3. Wavelet analysis 

Figures 8 and 9 are display both the original signal and its respective wavelet coefficients. 
Multilevel wavelet decomposition is carried out to divide the signal into two components: approximation 
coefficients and detailed coefficients. Single-level wavelet decomposition divides the signal into low-pass 
(approximation) components and high-pass (detailed) components. After the first level of decomposition, the 
approximation coefficients are obtained and further decomposed to result in second-level approximation 
coefficients (a2) and detailed coefficients (d2). This process is repeated for a total of 5 levels of 
decomposition, resulting in the coefficients d1, d2, d3, d4, d5, and a5. 

The wavelet analysis was carried out successfully and showed the waveform of the reactive power 
flow in the DFIG under healthy conditions. The waveforms in Figure 8 consist of various components of the 
reactive power absorbed by the DFIG, which validates its steady state. Next, the wavelet analysis was carried 
out with a fault introduced at 0.4 seconds. The waveform in Figure 9 shows the flow of reactive power in the 
DFIG with the fault. The various components of the reactive power consumed by the DFIG are shown with 
changes in amplitudes at 0.4 seconds, which validates the faulty state of the DFIG. To effectively assess the 
severity of the blade imbalance fault, it is expected that fluctuations in the amplitude component of the fault 
characteristic frequency of the stator current will be possible under the same wind speed conditions. 
Therefore, by comparing the variation in the frequency amplitude feature of I(t) with that in normal 
operation, it is possible to measure the severity of the blade imbalance defect. 
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Figure 8. Reactive power analysis of DFIG under normal condition 
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Figure 9. Reactive power analysis of DFIG under fault 


5. CONCLUSION 


mj 


2139 


This study analyses stator currents of a DFIG WT to identify blade imbalance conditions in the wind 
farm. The defect characterization methodology monitors the characteristic frequency Io™(t) to identify the 
severity of impeller imbalance flaws in real time, while wavelet analysis validates the results. The proposed 
approach analyses stator current characteristic extracts and diagnoses impeller imbalance failure in both 
steady-state and turbulent wind conditions. This method can accurately assess the severity of the blade 
imbalance fault by tracking variations in the stator current fault characteristic component's characteristic 


frequency amplitude under the same wind speed conditions. 
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